Pion Production by Protons on a Thin Beryllium Target 
at 6.4, 12.3, and 17.5 GeV/c Incident Proton Momenta 
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An analysis of inclusive pion production in proton-beryllium collisions at 6.4, 12.3, and 17.5 
GeV/c proton beam momentum has been performed. The data were taken by Experiment 910 at 
the Alternating Gradient Synchrotron at the Brookhaven National Laboratory. The differential -k^ 
and 7r~ production cross sections (d^a/dpdO,) are measured up to 400 mrad in 9^^ and up to 6 GeV/c 
in Ptt- The measured cross section is fit with a Sanford-Wang parameterization. 

PACS numbers: 13.85.Ni, 25.40.Ve 



I. INTRODUCTION 

A detailed understanding of the production of pions 
in proton interactions with nuclear targets is essential 
for determining the flux of neutrinos in accelerator based 
neutrino experiments. Flux predictions are particularly 
difficult for experiments using lower energy primary pro- 
ton beams such as MiniBooNE fit] (8 GeV) and K2K 
(12 GeV), where there currently exist large uncertain- 
ties in the pion production cross section data. Most of 
the existing data with proton beam energies in the 5 to 
20 GeV range were taken over 30 years ago using single 
arm spectrometers [1, 13, S S S IE] , but more recently the 
HARP Experiment at CERN has started to publish new 
data [§]. A global fit to the older data by Cho et al. [J| 
found a normalization discrepancy of ~15% between the 
various experiments and more recent fits [Tol . [Tl| have 
also found general inconsistencies. 

The high statistics data taken by Brookhaven experi- 
ment 910 (E910) provides an opportunity to revisit these 
old measurements with a modern, wide angle spectrome- 
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ter. The experiment covers a wide range of secondary mo- 
menta and angles with particle identification over most 
of this range. 

E910 has previously published tt^ production cross sec- 
tions for low momentum pions (0.1 to 1.2 GeV/c) on sev- 
eral different target materials including beryllium p^ . 
This analysis extends that earlier work to higher pion 
momenta. 



II. THE EXPERIMENT 

E910 ran for 14 weeks at the MPS facility in the Al 
secondary beam line of the BNL AGS in 1996. The to- 
tal momenta and directions of incoming beam protons 
were reconstructed using proportional chambers located 
upstream of the target. Data were taken at three beam 
momenta: 6.4, 12.3, and 17.5 GeV/c. Beam protons were 
identified by three Cerenkov counters along the beam 
line. 

A set of trigger counters (5*1, ST) and veto counters 
iyi, V2) located between the proportional chambers and 
the target were used to detect and constrain the trajec- 
tories of incoming beam particles. The trigger counters 
are shown in Figure [T] The trigger configuration for the 
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FIG. 1: The beam trigger counters and bullseye counter. The 
shaded area on the bullseye defines the beam veto region. A 
hit in this region indicates that the event is consistent with a 
non-interacting beam particle. 



data used in this analysis employs the "bullseye" counter, 
which was located 6.8 m downstream of the target. The 
bullseye counter consists of two pairs of scintillator slats, 
one pair of 14.6x30.5 cm slats placed along the verti- 
cal and the second of 40.6x7.6 cm slats aligned horizon- 
tally. Non-interacted beam particles consistent with the 
aperture defined by the veto counters arc entirely located 
within the intersection of one of the horizontal and one 
of the vertical slats as shown in Figure [H The minimum 
bias trigger for this analysis, known as the bullseye trig- 
ger, is defined to be the combination of the presence of 
beam (51 • VI ■ ST ■ V2) and the absence of a hit in the 
relevant slats of the bullseye counter. 

During the run, a variety of target materials were used. 
This analysis focuses on the proton-beryllium (Be) inter- 
action data sets. In the rest of the paper, only Be target 
data will be discussed. The Be target had a geometric 
cross section of 7.62x2.54 cm^ and was 1.84±0.04 cm 
long (~ 4.5% of an interaction length). The beam spot 
on the target was defined by the last veto scintillator, 
which had a 2x1 cm^ slot with semicircular ends. 

Reaction products from proton-beryllium interactions 
were measured with the spectrometer layout shown in 
Figure [21 The target was located 10 cm upstream of 
the time projection chamber (TPC) active volume. The 
EOS TPC [lil is 1.54 m long and is read out through a 
120x 128 cathode pad array. It ran with PIO gas at atmo- 
spheric pressure and a vertical electric field of 120 V/cm. 
The TPC was located in the center of the MPS mag- 
net, which had a nominal central field of 0.5 T along 
the vertical axis. Downstream of the TPC, inside the 
magnet, charged particle tracking was provided by three 
drift chambers (DCl-3). Each drift chamber consisted 
of seven wire planes: three x views (one staggered), two 
y views (staggered), and two views rotated from verti- 
cal by ±60°. A segmented threshold Cerenkov counter, 
with an aperture of 139.7x190.5 cm^, was located 4.8 m 
downstream of the target. The counter used 96 separate 
mirrors, a central 8x8 grid of small mirrors surrounded 
by 32 mirrors with a factor of four larger aperture, to 
direct Cerenkov light produced by particles traversing 
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FIG. 2: The E910 spectrometer layout. 



the nominally 1 m length of the counter to the same 
number of photomultiplier tubes located at the top of 
the counter. The radiating medium was Freon 114. The 
time-of-flight (TOF) wall was located 8 m downstream of 
the target. It consisted of 32 scintillating counter slats, 
each 15.2x178 cm^, arrayed in a flat panel of approxi- 
mately 488x178 cm^. The typical TOF resolution was 
^160 ps [l3|- Two more drift chambers (DC4-5) were 
located beyond the TOF wall at 9.6 and 10.1 m from the 
target. 



III. DATA ANALYSIS 

The differential cross section for pion production as a 
function of pion momentum, p, and zenith angle, 0, is 
given in the lab frame by 

dpdn^^' Na^P e a{p,e) Ap2ttAcos9 Nbeam 

where A is the target mass number (Abe — 9.01 g/mole), 
Na is Avogadro's number, £ is the target length, p is 
the target density [pse — 1.848 g/cm"^), e is the trigger 
efficiency, a(p, 9) is the geometrical acceptance and cut 
efficiency measured (with Monte Carlo simulation) as a 
function of p and 9, Ap27rAcos0 is the area of the bin 
in {p ri) space, iV^ (p, 9) is the number of pions observed 
in the bin, and Nbeam is the total number of protons 
incident on the target. 

The data are binned in 400 MeV/c wide bins in mo- 
mentum. In the 12.3 and 17.5 GeV/c data sets, the 9 
bin width is 60 mrad from to 360 mrad, and for the 
6.4 GeV/c data set, the 9 binning is 100 mrad from to 
400 mrad to accommodate lower statistics. 



A. Event Selection 

For this analysis, all events - both those used to de- 
termine Nbeam and A'^ - must satisfy the following cri- 
teria: the beam particle must be successfully tracked in 
the beam proportional chambers; the beam Cerenkov re- 
sponse must be consistent with the proton hypothesis; 
the reconstructed primary vertex must be consistent with 
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the target z position and lie within the open apertur of 
the V2 counter. 

A completely unbiased beam trigger is used to deter- 
mine the the number of beam protons, Nheam- This trig- 
ger fired on the presence of beam only, without requiring 
an interaction. The trigger was heavily prescaled; by fac- 
tors of 32 to 64 depending on the run. The number of 
beam protons is determined by counting all beam trigger 
events and multiplying by the prescale factor. As a cross 
check on the prescale, the prescale factor is computed in 
the data by taking the ratio of buUseye trigger events 
which also have a prescaled beam trigger to all bullseye 
triggers. In all data sets this number agrees with the 
input prescale to well within the statistical error on the 
ratio. 

Candidate events are required to have a bullseye trig- 
ger, and a reconstructed vertex position consistant with 
the reconstructed beam track and within the volume of 
the target. The reconstructed secondaries must be con- 
sistent with a single inelastic collision (i.e. not consistent 
with an elastically scattered beam particle and the sum 
of the secondary particle momenta must not be greater 
than the beam momentum). 



B. Track Selection 

All candidate tracks must point back to the primary 
interaction vertex, they must have a minimum of 20 hits 
in the TPC, and the TPC dE/dx calculation must have 
been successful. 



C. Trigger Efficiency 

The bullseye trigger efficiency, e, was determined by 
using a sample of totally unbiased beam trigger events. 
The denominator of the efficiency is the number of beam 
trigger events with at least one secondary track, and the 
numerator is the subset of those events that also have a 
bullseye trigger. The measured trigger efficiency for each 
beam momentum data set is listed in Table HI The error 
on the trigger efhciency is due only to statistics. 

In addition to the flat trigger inefficiency, which is 
largely due to other tracks in the event that pass through 
the bullseye, we studied the possible inefflciency as a 
function of pion kinematics {i.e. for a particular angle 
and momentum the track in question may have a non- 
zero probability of passing through the bullseye). This 
effect was determined to be small (<2%) in all bins and 
mostly affects the lowest angle bin. The bin-to-bin inef- 
ficiency is accounted for as a systematic on a bin-by-bin 
basis. 



Data Set 


Beam Protons 


Trigger Efficiency 


6.4 GeV/c 


93,632 


100.il. 1% 


12.3 GeV/c 


745,216 


96.8±0.6% 


17.5 GeV/c 


2,576,352 


89.6±0.6% 



TABLE I: Number of protons on target and trigger efficiency 
for each data set. 

D. Particle Identification 

Three particle identification (FID) systems are used 
to distinguish pions from other secondaries in this analy- 
sis: the TPC, timc-of-flight wall, and threshold Cerenkov 
counter. The TPC uses energy loss (dE/dx) to distin- 
guish between different particles types. The pion thresh- 
old in Cerenkov counter is about 2.5 GeV/c while the 
kaon and proton thresholds are both outside the momen- 
tum range of this analysis (9.0 and 17.1 respectively). 
The TOF wall performance is discussed in ref . 

For each of the PID detectors, a residual is formed be- 
tween the true response and the expected response for the 
different particle hypotheses (e, tt, K, p). By construc- 
tion, the residual for the correct hypothesis is Gaussian 
distributed with a mean of zero and a width of one. 

Pion candidate tracks are divided into three groups as 
a function of momentum. The groups are defined by the 
capabilities of the three PID systems. In each momen- 
tum region the pion residual for the primary PID system 
is plotted for all candidate track in each {p, 9) bin. These 
residual distributions are then fit to a unit Gaussian to 
determine the pion yields. Information from the other 
two PID systems is used as a discrete cut (described be- 
low) where applicable. 

In the range of 0.4 to 1.2 GeV/c the TPC dE/dx is 
used as the primary PID system (the dE/dx distributions 
for different particle types can be found in Figure 3 of 
Ref. [13] )■ The pion-proton and pion-electron separation 
is excellent over this entire range. The pion-kaon residu- 
als start to overlap at the 3 cr level by about 0.8 GeV/c. 
Kaon production is smaller than pion production by a 
factor of ~20, and the residual overlap is not complete. 
Therefore, kaon contamination represents an error on the 
pion yield of much less than 5%. The assertion that the 
kaon contamination is small is affirmed by the quality of 
the fits to the pion residuals. The TOF PID system has 
a small geometrical acceptance in this range, but if TOF 
information is available for a track, the TOF pion resid- 
ual is required to be within ±3 a of zero. The Cerenkov 
system has no ability to separate pions from kaons and 
protons in this region and therefore is not used. 

Above 1.2 GeV/c, the primary PID system is the TOF. 
In the TOF system pion-kaon separation is good (> 3ct) 
up to ~3 GeV/c and pion-proton separation is good up 
to ~5.4 GeV/c. Above 2.8 GeV/c all pions should pro- 
duce a robust Cerenkov signal, while kaons and protons 
do not. In this region, the Cerenkov pion residual must 
be within 3 cr of zero for all pion candidate tracks. Both 
below and above 2.8 GeV/c, there is a requirement that 
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TPC dE/dx pion residual be within 3 cr of zero. Below 
2.8 GeV/c, this cut is useful for separating pions from 
electrons. Above ^^3 GeV/c, the relativistic rise of the 
pion in dE/dx provides additional pion-proton separa- 
tion. 



E. Acceptance and Efficiency 

A large sample of Monte Carlo events (~680,000) was 
used to determine the pion geometric acceptance and cut 
efficiency, by taking the ratio of generated to accepted 
TT"*" and 7r~ in each (p, 6) bin. The product of acceptance 
and reconstruction efficiency was determined by binning 
in the generated particle momenta. The effects of finite 
momentum and angular resolution on the reconstructed 
spectra were separately evaluated as systematic errors. 
On a bin-by-bin basis, the size of the smearing uncer- 
tainty is generally less than 5%. 



IV. RESULTS 

The inclusive pion production cross section in proton- 
beryllium interactions is calculated from candidate pion 
tracks, binned in p and 9. Each momentum bin spans 
0.4 GeV/c. The first bin begins at 0.4 GeV/c and the 
final bin ends at 5.6 GeV/c. Six zenith angle bins, rang- 
ing from to 360 mrad, are used with the 12.3 and 
17.5 GeV data sets. For the much smaller 6.4 GeV data 
set only four angular bins are used covering the range 
to 400 mrad. The angular bin spacing is uniform in 6, 
and the bin centers are reported as the average cos 61 of 
the bin. Tables IIIIl IIVI and |V] list the measured value of 
the TT"*" and 7r~ cross sections in each (p, 9) bin, in units 
of mbarns/[(GeV/c) steradian]. 

The errors reported in Tables [ml [IV] and |V] include 
contributions from data statistics, Monte Carlo statistics 
and bin-by-bin cross checks. These checks include: PID 
studies, bin migration studies, and a study of buUseye 
trigger inefficiency as a function of p and 9. On average, 
the largest systematic contribution comes from a com- 
parison of the cross section of tt^ production measured 
with and without the PID cuts. To first approximation, 
all negative secondary tracks are pions. Therefore, it is 
possible to measure the ■k~ cross section without using 
PID cuts. This PID-free cross section should agree well 
with the full PID based analysis, and any areas where 
the PID cut efficiency, as calculated in the Monte Carlo, 
is not a perfect match to the data would be highlighted. 
One expects that deviations might appear in the tran- 
sition regions such as just above 1.2 GeV/c where the 
primary PID system switches between TPC and TOP, 
and around the pion threshold in the Cerenkov counter 
at ~2.8 GeV/c. Bin-for-bin the difference between the 
two analyses is taken as a measure of PID systematic 
error which is applied to both the 7r~ and tt"*" error anal- 
ysis. The PID systematic is typically less 5%, but can be 



larger in the PID overlap regions. 

The error in each bin (reported in Tables lllli IIVI and 
W\ is the quadratic sum of all systematic contributions 
and the statistical error. In addition, there is an over- 
all normalization uncertainty which should be applied 
equally to all bins. This error is due in part to a 2% 
uncertainty on the measurement of the target thickness, 
and to the uncertainty in the trigger efficiency shown in 
Table m The total normalization error is estimated to be 
less than 5%. 



V. FIT TO THE DATA 

Modeling the pion production cross section as a func- 
tion of beam momentum, secondary particle momen- 
tum, and secondary particle angle is of interest for in- 
put to Monte Carlo simulations. For this purpose, we fit 
the results of the previous section with a Sanford-Wang 
function [Tsj. which describes the inclusive double differ- 
ential pion production cross section in proton-beryllium 
interactions. 

Sanford-Wang fits have been used in recent years to 
describe the inclusive pion production cross section for 
low energy neutrino experiments (lo| . [ll[, and in the 
more distant past for global fits to low energy inclusive 
pion production data 5] . The functional form of the pa- 
rameterization was developed empirically, based on data 
with incident proton momenta between 10 and 70 GeV/c, 
therefore fitting the 6.4 GeV/c data provides a useful test 
of its range of validity. The explicit form of the Sanford- 
Wang parameterization used in this analysis is 



SW = 



exp 



Pb - I 



^cf"" - C6 9^ {p„-cr Pb cos'^* 9^ ) 



(1) 



Pb 



where is the momentum of the pion (in GeV/c), pf, 
is the momentum of the beam proton (in GeV/c), is 
the production angle (in radians) of the pion in the lab 
frame, and the Cj are parameters to be obtained by a fit 
to the data. 

To fit for the parameters, Ci through Cs, we use the 
following function: 



E 



[Nj X SW 



\ dpcin J . 



i 



(2) 

where i spans all data points in p^r and , and j spans the 
three beam momentum settings; Nj is a normalization 
term for the jth beam momentum; SWi is the Sanford- 
Wang function (Eqn [T]) averaged over and 9^^ in the 
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ith bin: 

SWi = - — -4 7- X (3) 

/ / 5M^(pb,p^,6'^)sin6'^(ip^d6'^; 

(cPu / dp(Kl)i is the measured cross section in the ith bin, 
and (Ti is the measurement error on bin i including both 
systematic and statistical errors (the normalization error 
is not included) . The normalization uncertainty of each 
incident proton momentum data set is handled with the 
terms N , which add to the relative to the normaliza- 
tion error, ctjv- For all incident proton momentum data 
sets, ctat is taken to be 5%. 

The best fit parameters for the Sanford-Wang fit to 
the combined 6.4, 12.3, and 17.5 GeV/c 7r+ (tt^) data 
are shown in Table |TT1 The per degree of freedom for 
the 7r+ (tt") fit is 1.93 (1.59). For both fits, the values of 
the normalization terms at the best fit point are within 



1 cr of zero: 1.03 (1.01), 1.01 (1.01), and 0.97 (0.98), for 
the 7r+ (tt") fit in order of increasing beam momentum. 
The parameters are highly correlated, so their errors are 
given as a full covariance matrix, which for the 7r+ fit is: 



Parameter 


7r+ Data Fit 


TT Data Fit 


cl 


258.2 


249.3 


c2 


1.018 


1.066 


c3 


2.953 


3.311 


c4 


2.204 


1.188 


c5 


1.782 


1.017 


c6 


5.136 


5.127 


c7 


7.706 X 10"^ 


6.459 X 10"^ 


c8 


14.64 


10.22 


x' 


323 


268 


Na.o.f. 


167 


169 



TABLE II: Best-fit parameters for Sanford-Wang fits to the 
combined 6.4, 12.3, and 17.5 GeV/c production data sets in 
■k'^ and 7r~. 



/ 210. 


0.587 




10.6 


5.69 


1.37 




6.36 


0.124 




57.0 


0.587 


1.67x10- 


3 


2.98x10-2 


1.60x10-2 


3.85x10 


-3 


1.79x10-2 


3.48x10" 


-4 


0.160 


10.6 


2.98x10" 


2 


0.541 


0.289 


6.97x10 


-2 


0.323 


6.30x10- 


-3 


2.89 


5.69 


1.60x10- 


2 


0.289 


0.156 


3.74x10 


-2 


0.173 


3.38x10- 


-3 


1.55 


1.37 


3.85x10- 


3 


6.97x10-2 


3.74x10-2 


9.04x10 


-3 


4.17x10-2 


8.14x10- 


-4 


0.374 


6.36 


1.79x10- 


2 


0.323 


0.173 


4.17x10 


-2 


0.195 


3.77x10- 


-3 


1.73 


0.124 


3.48x10- 


4 


6.30x10-3 


3.38x10-3 


8.14x10 


-4 


3.77x10-3 


7.40x10- 


-5 


3.39x10 


\ 57.0 


0.160 




2.89 


1.55 


0.374 




1.73 


3.39x10- 


-2 


15.7 



and for the tt fit the covariance matrix is: 



/ 197. 


0.912 




8.21 




2.03 




1.40 




6.33 




7.68x10- 


2 


38.3 


0.912 


4.56x10 


-3 


3.82x10- 


2 


9.50x10- 


3 


6.49x10- 


-3 


2.97x10- 


2 


3.58x10- 


4 


0.178 


8.21 


3.82x10 


~2 


0.345 




8.52x10- 


2 


5.85x10- 


-2 


0.265 




3.22x10- 


3 


1.60 


2.03 


9.50x10 


~3 


8.52x10- 


2 


2.14x10- 


2 


1.45x10- 


-2 


6.56x10- 


2 


7.98x10- 


4 


0.398 


1.40 


6.49x10 


-3 


5.85x10- 


2 


1.45x10- 


2 


1.00x10- 


-2 


4.52x10- 


2 


5.48x10- 


4 


0.273 


6.33 


2.97x10 


-2 


0.265 




6.56x10- 


2 


4.52x10- 


-2 


0.208 




2.49x10- 


3 


1.24 


7.68x10-2 


3.58x10 


-4 


3.22x10- 


3 


7.98x10- 


4 


5.48x10- 


-4 


2.49x10- 


3 


3.04x10- 


5 


1.50x10 


V 38.3 


0.178 




1.60 




0.398 




0.273 




1.24 




1.50x10- 


2 


7.54 



Figures H [51 and [7] show the 6.4, 12.3, and 17.5 GeV/c 
proton beam momentum tt"*" data sets compared with 
the fit result overlaid. The corresponding distributions 
for TT~ data are shown in Figures IH [6l and [51 

Fitting the data to a parameterization like Sanford- 
Wang allows for the calculation of the pion production 
cross section as a function of incident proton momentum 
and for a wide range of pion angles and momenta. In par- 
ticular, this fit can be used to generate the primary pion 
production for experiments such as MiniBooNE with a 
primary beam momentum (8.9 GeV/c) that was not di- 



rectly studied, and the fit covariance matrix can be used 
to calculate the uncertainty in that primary production 
model. 

The fit is in reasonably good agreement with the data 
in all but the lowest angular bin (42 mrad). In this region 
the fit appears to be systematically above the data, espe- 
cially in the case of 17.5 GeV/c protons and 12.3 GeV/c 
protons with secondary tt momenta in the 1 to 3 GeV/c 
range. This may be due a deficiency in the Sanford-Wang 
parameterization, or it could be from a systematic and 
unaccounted for inefficiency in the low angle region. 
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e (mrad) p (GeV/c^) 7r+ 7r~ 6* (mrad) p (GeV/c^) 
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1.4 


131.6 


35.2 


87.5 


20.4 




1.4 


100.1 


26.8 


46.5 


11.7 


1.8 


131.4 


27.7 


42.4 


15.5 




1.8 


79.3 


21.2 


17.9 


10.0 


2.2 


95.3 


22.7 


43.8 


14.9 




2.2 


76.7 


18.4 


21.0 


7.4 


2.6 


59.0 


85.2 


14.7 


84.0 




2.6 


19.6 


19.6 


13.4 


18.2 


3.0 


48.5 


18.3 


48.6 


14.4 




3.0 


13.2 


6.8 


7.9 


3.2 


3.4 


27.6 


27.7 


11.2 


25.7 




3.4 


12.1 


6.7 


17.9 


5.0 


3.8 


5.6 


7.2 


5.7 


8.1 




3.8 


— 


— 


6.6 


2.7 


4.2 


6.0 


6.1 


6.0 


4.1 














0.6 


111.6 


9.9 


68.7 


7.6 


353 


0.6 


84.9 


8.3 


65.1 


6.9 


1.0 


109.9 


10.3 


53.5 


7.2 




1.0 


56.9 


6.5 


40.5 


5.4 


1.4 


76.8 


22.9 


26.0 


9.3 




1.4 


44.5 


44.6 


20.7 


41.1 


1.8 
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TABLE III: Pion production cross sections for 6.4 GeV/c protons on Be. 




FIG. 3: Inclusive p-Be n production cross section data and fits 
vs. 7v^ momentum, at 6.4 GeV/c incident proton momentum. 
Fits are defined in Table HT] 



FIG. 4: Inclusive p-Be production cross section data and fits 
vs. Tv" momentum, at 6.4 GeV/c incident proton momentum. 
Fits are defined in Table HH 
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TABLE V: Pion production cross sections for 17.5 GeV/c protons on Be. 
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FIG. 5: Inclusive p-Be n production cross section data and fits 
vs. tt"*" momentum, at 12.3 GeV/c incident proton momentum. 
Fits are defined in Table HH 



FIG. 6: Inclusive p-Be tt" production cross section data and fits 
vs. 7r~ momentum, at 12.3 GeV/c incident proton momentum. 
Fits are defined in Table HH 
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FIG. 7: Inclusive p-Be 7r+ production cross section data and fits 
vs. tt"*" momentum, at 17.5 GeV/c incident proton momentum. 
Fits are defined in Table HI 



FIG. 8: Inclusive p-Be tt" production cross section data and fits 
vs. 7r~ momentum, at 17.5 GeV/c incident proton momentum. 
Fits are defined in Table HIl 



